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Abstract — This paper presents identification of a non-linear 
system for robot flexible arm as a black-box using Volterra 
series model. In most of the methods using Volterra series, too 
many number of parameters are estimated which is in contrast 
with the simplification principle in system identification. In this 
paper, for system identification of robot arm, kernels of Volterra 
series are developed using orthogonal functions that would 
reduce the number of parameters and consequently decreasing 
the identification time and increasing the accuracy and speed of 
convergence. At last, results of system identification of flexible 
robot arm based on expansion of Volterra series kernels through 
kautz orthogonal function are investigated with the results of the 
other method in literatures. 

Index Terms — Volterra Series, kautz Orthogonal Functions, 
non-linear System Identification, Optimization, flexible robot 
arm 

I. INTRODUCTION 

Despite numerous methods for identification of linear 
systems, there are a few methods for non-linear systems 
(Dalirrooy Fard and Karrari ,2005). One of the most 
comprehensive methods for modelling non-linear systems is 
Volterra Series. Generally speaking, as modeling with 
Volterra series would not interfere with the hidden state 
variables of the system, the possibility of specifying the 
dynamic behavior of the system is provided without 
measuring the state variables (Minu and Jessy John ,2012). 
Volterra Series produces a model for non-linear behaviour 
of the system like Taylor series, but in Taylor series output 
is highly dependent to the input in a certain instant Whereas 
In Volterra series, output is dependent to the input in all the 
instants (S .Parker et al.,2001). The main advantage of the 
Volterra series lies in its totality, so many non-linear systems 
can be modeled by it (Mahmoodi, et al.,2007). This series 
has high number of parameters; therefore it is named a non- 
parametric model. Series and theory of Volterra were first 
introduced by Vito Volterra in 1887 (Volterra ,1887; 
Volterra ,1958). After that, this series was widely used in 
identifying the non-linear systems, where a sample 
application of this series is presented in (Alper ,1965). 
Literatures about this issue are numerous and most of them 
try to identify a certain system in a certain application or 
presenting a new method for calculating the kernels of this 
series from which ref (Meenavathi and Raj esh, 2007) is a 
good example that uses Volterra series to remove Gaussian 
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noise. In (Dalirrooy Fard and Karrari ,2005), Volterra series 
is used for modeling and identification of the synchronous 
generator. Ref (George-Othon et al., 1999) presents the 
specification of Volterra series kernel with multiple 
methods. Ref (Dodd and Harrison, 2002 ; Yufeng Wan et 
al., 2003) use Volterra series for estimation of high rank 
kernels through reproduced Hilbert space. 

In this paper in section 2, Volterra series is described and in 
section 3, the method for calculation of kernels or factors of 
Volterra series are presented. Laguerre orthogonal function is 
described in section 3-1 and expansion of kernel with 
orthogonal function is presented in section 3-2. Flexible robot 
arm is illustrated in section 4 and section 5 presents the results 
of simulation. 

II. VOLTERRA SERIES 

Volterra series is generalized idea on convolution in linear 
system. The relation between input and output in a linear and 
memory less system is expressed as follow: 
y(0=h.u(t) (1) 

in which output is only dependent to the input in the instant 
t. In a linear, casual, time-invariant, with memory the 
relation between input and output is as follow (Doyle et 
al.,2001): 

oo 

y\( n ) = y %l,K( m i )u(n-m i ) (2) 

i= 0 

The relation between input and output in a non-linear and 
memory less system of rank two is expressed as follow: 

y 2 (t) = h 2 .u\ty (3) 

The relation between input and output in a linear and with 
memory system of rank two is expressed as follow: 
y 2 (n) = u(n).u(n) + u(n).u(n - 1) + ... 

+ u(n).u( 0) + u(n - 1 ).u(n - 1) + ^ 

u{n - 1 ).u(n - 2) + ... + u(n - l).u(0) + .... 

+ u(0).u(0) 

Where it can be written 

00 00 

y 2 (n) = ^ u(n)u{n - j ) u(n - 1 )u(n- i ) + ... 

j = o ;= o 

(5) 

00 00 

= -!>(«- 7 ) 

7=0 i= 0 

Now we can add a proper weight to the series of equation (5) 
without disturbing the whole problem: 
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oo oo 

V V i 


y 2 (n) = XX h 2 (*> /)•“(” - 0-w(ra - 7) (6) 

J =0 z '=0 

For Volterra series of non-linear system with rank n: 

y - % + y\ + y 2 +— + y n (7) 

where y- x is calculated like equation (2) and (5) and y 0 is a 
constant value (e.g. signal average). Now, substituting y* in 
(7) gives: 


y(n) = y 0 + Y^h x (,m x ) u (n-m x ) 


+ 


m l =0 


( 8 ) 


y y /? 2 (m l ,m 2 )u(n-m 1 )u(n-m 2 ) + ... 

m 1 =0ra 2 =0 

h n is rank n kernel of Volterra. In fact hi(t) shows impulse 
function of linear system and h n (t) is impulse function of a n- 
dimensional non-linear system. Equation (8) has too many 
unknowns that calculation of it is time-consuming. For 
solving this problem we have to consider that there may not 
be link between output and former inputs. Besides, series 
can be shown with kernels of rank 1 to L by approximation. 
From these points, Volterra series is changed as follow: 

M - 1 

y(n) = y 0 + Yj h \ (rn l )u(n-m l ) + 

nil =0 

(9) 

M-l M - 1 

V V h 2 (m, ,m 2 )M(n- m ] )u(n- m 2 ) 

m l =0 m 2 =0 

M (memory of system) and L (rank of system) are almost 
chosen by trial and error method. 


III. . CALCULATING THE KERNELS OF VOLTERRA 
SERIES 

Calculation of Volterra series factors is the main issue in 
solving the equation and identifying the systems using this 
series. Identifying non-linear systems by Volterra series is in 
fact identifying the kernel of such series. There are many 
methods for this purpose such as recursive algorithm (Doyle 
et al.,2001), Gaussian input (Doyle et al.,2001), gradian- 
based search (Suleiman and Monin ,2008) , cross-correlation 
method (Lee and Schetzen ,1965), method of reproducing 
kernel Hilbert space (Dodd and Harrison, 2002 ; Yufeng 
Wan et al., 2003) and etc. In addition, the factors of Volterra 
series can be obtained directly by expanding the orthogonal 
functions that leads to less complexity and more accuracy. 
The main importance of such method is that if we rewrite 
kernels of Volterra series as an expansion of orthogonal 
function, series is changed into a linear regression equation 
(Moodi ,2008) In this paper, number of kernels is also 
considered symmetrical: 

h s „ ym (k v k 2 ,...,k n ) = ± ( 10 ) 

Due to complexity of issues for identifying the kernel of 
Volterra series and writing limitation, we only illustrate the 
method for identifying kernel by Laguerre orthogonal 
function. 

IV. Kautz ORTHOGONAL FUNCTIONS 
Orthogonal Katz function 

The two- parameter Katz function is defined as follow: 


p GO 



'-cz 2 -b £j(c - l) + z -b 1 

z 2 -E & (c — l)z - c 

z 2 -b b[c - l3 z - c 


(ID 


f ^ Cl - c : )ztz - b) 

02n - lCzJ = z’ + bCc - tfz - c 
'-cz 2 -f b (c - 1) -f z -b r n 1 

z 2 -\-b(c — l)z — c 


(123 


In which a, b are real value constants related to the poles of 
Katz function. 

b = (fi +/JX1 (133 

c = -Pfi (143 

Simultaneous and optimal choice of c and b is being 
investigated. So, to obtain the best choice, these parameters 
considered constant by applying a certain criteria. 


V. EXPANSION OF KERNELS WITH ORTHOGONAL 
FUNCTIONS 

As mentioned, expansion of orthogonal basic functions is an 
estimation method that reduces the number of parameters 
and consequently reduces the complexity and increases the 
accuracy and also changes the series to a linear regression. 
Now factors are calculated as below(Moodi and Bustan 
, 2010 ): 

p 

h Mi ) = Tj a iAS m i) 

*i=i 

(17) 

P P 

h 2 (m t , m 2 ) = X X a ¥2 A ( m i M 2 ( m 2 ) 

*i =1 h =1 

Where ® k (n) is the kth orthogonal function at point n. so, the 
relation between input and output in a Volterra series of rank 
two is as follow: 

y («) = y 0 + L a,, (n) + L X a , h ^ A ^ n8 > 

ly — 1 ly —1 = 1 

In this equation factors are symmetrical ( a ^ = a ), 
therefore: 

M - 1 

S k (n)=Y,<t>k(m l )u(n-m l ) (19) 

m, =0 

First, M and L and P are chosen. Although the actual values 
of these parameters are too high, if values with lesser degree 
have an agreeable error, there is no need to increase the 
values. Finally the equation (18) is written as: 
y = U0 + e (20) 

where y is output vector, 0 vector of factors and e noise or 
model error, U and known matrix of S values. Now 
different identification methods are used in a way that 

modeling error || y — y|| 2 is minimized. One of such methods 

is linear least square that best parameters are calculated as 
follow: 

P L s ~ (U T U)~ 1 U T y (21) 
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The above method has some problems and that is inversing 
from U T U. For example maybe det(U T U) is not zero but be 
near zero. In this case, matrix is called improper and the 
reason is the dependency between rows and columns of U T U 
that identification would be with less accuracy and high 
error. There are some reasons for this (kararri,2009): 

I. Input is not exciting enough and not able to excite 
all the modes of the system. 

II. The considered rank of model is more than actual 
rank of system. 

III. System is identified in closed-loop operation. 

For compensating this problem in least square method, 
inversing U T U is not allowed and other techniques are 
employed. In this paper SVD is chosen for this purpose. 
Besides, equation (22) is a criterion for feasibility 
assessment of model under the name of forecasted error 
(Barve ,2004): 


e 


p 



}’;(]) ~ YU)) 2 


7=1 


N - 1 


Zuu )) 2 

7=0 


( 22 ) 


where P is number of outputs, N number of used data and 
y.(y) actual output of i at instant j and y t (j) output of the 


model of i at instant j . 


expansion coefficients of the model structure is obtained in 
IS way. 



k 


Figure 2. Output waveform of actual system and detected 
output 

To validate and display the success of identification made, 
the criterion of Mean Squared Error (MSE) and the 
Normalized Root Mean Square error (NRMS) are calculated 
as follow: 

1 N 

MSE = — Y,e(kf = 5.0891 x 10 04 
N k = 1 


hr( y (/c)-mf 


NRMS 


= 0.0820 


2>w 2 
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VI. FLEXIBLE ROBOT ARM 

Robot arm is a system with high applicability in industrial 
systems which are shown in Figure 1 (Ferreira and Serra 
,20 12). Flexibility of the robot arm is in two types: 

• Flexibility due to non-linear variations in Drive 
system. 

• Flexibility due to straining the robot arm. 

Data of input and output of the system are categorized into 
two groups, one group for training that the model is 
extracted with and the other one for experimenting the 
model. 



VII. SIMULATION RESULTS 
1 . Simulation result 

To identify the flexible robot arm to the structure of Volterra 
series model developed by Katz function will be considered 
according to Equation 16. First, we need to do optimal 
estimation of Katz orthogonal function based on the real 
measured input and output available in the paper of ALEX. 
According to the model structure, input and output data, the 


VIII. CONCLUSION 

In this paper, a comprehensive model based on Volterra 
series was presented for identification of non-linear system 
of flexible robot arm. For reducing the complexity and 
increasing the accuracy of the model, kernels of Volterra 
series were rewritten as an expansion of orthogonal 
functions and optimal values of orthogonal functions pole 
and number of expansion. For preventing inversing 
problems and impropriety of regression matrix, were 
compared and it was shown that identification with 
developed Volterra series by kautz function has less error 
and more accuracy comparing to other method. In addition, 
there are other methods of optimization for better adjustment 
of orthogonal functions which can improve the identification 
process. This is hoped to be done in future works. 
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